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Abstract 

We consider the prompt hadroproduction of J/if), ip' and the T(IS,2S,3S) states 
caused by the fusion of a symmetric colour-octet state, (gg)ss, and an additional gluon. 
The cross sections are calculated in leading-order perturbative QCD. We find a consid- 
erable enhancement in comparison with previous perturbative QCD predictions. Indeed, 
the resulting cross sections are found to be consistent with the values measured at the 
Tevatron and RHIC, without the need to invoke non-perturbative 'colour-octet' type of 
contributions. 



1 Introduction 

It is not easy to describe the hadroproduction of J /if) mesons within a perturbative QCD 
framework. The problem is that, due to the J p = 1~ quantum numbers, it is not possible to 
directly form the colourless J /if) meson by gluon-gluon fusion. The simplest possibility is to 
produce a colour-octet quark-antiquark pair (gg — > cc) and then to emit an additional gluon, 
which carries away the colour, as shown in Fig.^i. This is often referred to as the colour- singlet 
mechanism (CSM). However the corresponding cross section is suppressed by the small QCD 
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Figure 1: (a) The 'bleaching' gluon subprocess used in the original 'colour- singlet' perturbative 
QCD estimates of prompt J /if) hadroproduction. (b) The perturbative QCD mechanism studied 
in this paper. In each case the subprocess gg — > J/if> g is shown in bold. 



coupling a s , and by the additional phase space factor associated with the extra gluon emission. 
As a result the LO QCD prediction pQ is found to be about an order of magnitude lower than 
the experimental yield of J/if> mesons. 

An alternative and more phenomenological approach is provided by the Colour Evaporation 
Model (CEM) [2|. Here the quarkonium production cross section is an (a priori unknown) 
fraction of the QQ heavy quark cross section integrated over the rnqq invariant mass up to the 
threshold for producing a pair of the lightest heavy flavour mesons. There are no constraints 
on the colour or spin of the QQ pair, the transition from colour-octet QQ to colour- singlet 
quarkonium is assumed to take place by the 'evaporation' of soft gluons. The fraction of QQ 
pairs that materialise as a particular quarkonium state, f c5 ^j/^ for example, is assumed to be 
universal and is adjusted to give the best fit to existing data. 1 Despite its phenomenological 
success, the CEM has no firm theoretical foundation. In practice one would not expect the 
evaporation of soft gluons to take place independently of the particular collision environment, 
and there is no reason why such soft interactions would not modify the quarkonium production 
properties, and in particular its collision energy dependence. This is precisely what happens in 
the theoretically more rigorous formalism for quarkonium production that we introduce below. 

The currently most popular and widely used description of quarkonium production is based 
on a nonrelativistic QCD (NRQCD) effective field theory approach [3], which retains features 
of both the CSM and CEM. Here QQ pairs are produced via a hard partonic (short-distance) 
subprocess, in both colour- singlet and colour-octet states, and non-perturbative universal ma- 
trix elements describe the (large-distance) transition of the QQ pair into particular quarkonium 

practice, it will of course depend on the parameters and pdfs used in the calculation of the QQ cross 
section. 
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states. 2 . The cross section of quarkonium H production is written in the schematic form 



da{H) 



Y,da{QQ[n\) (0 H [n}), 



(1) 



n 



where n denotes the set of colour and angular momentum quantum numbers of the QQ pair, 
and a is the cross section of the QQ pair production in a hard subprocess. The non-perturbative 
transition from the QQ state n into the quarkonium state H is described by a long-distance 
matrix element (0 H [n\). These matrix elements are taken as parameters to be determined by 
fits to experimental data. In this way, it is possible to compensate the low value of the hard 
subprocess cross section, da{QQ[n]), by the large fitted value of the matrix element {0 H [n}). 
Despite some phenomenological success, a detailed proof of the factorisation formula (HJ) is 
lacking, and in particular it is expected to break down at small values of the quarkonium 
transverse momentum. Predictions for the total quarkonium cross section in the NRQCD 
approach must therefore be treated with caution. 

Since the mass of the J/ ip meson is not particularly small, it would be desirable to be able 
to describe the QQ — > H transition within a perturbative QCD framework. In other words, we 
would like to consider explicitly an extra gluon exchange (which was hidden in the value of the 
non-perturbative matrix element (0 H [n}) in the case of the colour-octet mechanism, i.e. where 
the J/ip meson is formed from the colour-octet QQ pair after some non-perturbative interaction 
described by the (0 H [n\) matrix element. The corresponding lowest order in a s diagrams are 
shown in Fig. El In comparison with the gg — > QQ amplitude, the contributions of Fig. El 
contain an extra loop factor with a small coupling a s . Here we investigate the possibility that 
this suppression is compensated by the large number of graphs where the additional (third) 
gluon, needed to form the J/ip, is absorbed by different parton-spectators. As viewed from 
the collinear approximation, the amplitude shown in Fig. [TJd, and Fig. |2J corresponds to the 
next-to-next-to-leading order (NNLO) contribution to the cross section (as compared with LO 
99 ~ > Xc production). However the second t-channel gluon may be absorbed by any parton 
spectator; that is the amplitude is enhanced by the parton multiplicity, n oc logs. Therefore 
it may be considered as the LO amplitude in the BFKL approach. The aim of this paper 
is to evaluate the numerical size of this enhanced 'NNLO' contribution, and to see if it can 
remove the large discrepancy between the perturbative QCD prediction and the inelastic J/ip 
hadroproduction data. 

In Section 2 we calculate the amplitude of the process shown in Fig. El using the non- 
relativistic QQ — > J/tp vertex that was proposed in Ref. [7j. Due to the non-relativistic wave 
function of the J/ip, there is practically no integration over the quark loop in Figs.|2k,b. Indeed, 
the J /ip vertex (i.e. the J/ip wave function integrated over the relative momenta of the charm 
quarks), together with the two nearest c-quark propagators is 



where the index v corresponds to the J/ip polarization vector; and = Q^/2 and m are the 
4- momentum and the mass of the c-quark (Q M is the momentum of the J/ip)- The constant g 

2 See, for example, 0JE1E] and references therein, for a more detailed discussion of the situation, and of the 
history of the subject. 



gtft + m)^ 



(2) 
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Figure 2: Lowest-order perturbative QCD diagrams for J/ijj hadroproduction via gluon-gluon 
fusion with an additional gluon. 



may be expressed in terms of the electronic width T J ee of the J/ip — > e + e decay 

9 ~ 64^' ^ 

where M = 2m is the mass of the J/ip meson, and the electromagnetic coupling a = 1/137. 

In Section 3 we compute the total cross section of prompt inelastic J/ip hadroproduction 
at collider energies. That is, in that section, we neglect the additional J/ip yield coming from 
6-quark or Xc decays. We find that the LO result agrees with the Tevatron data rather well. 
The transverse momentum distributions of prompt J/ip, ip' and the T(1S,2S,3S) production 
are presented in Sections 4 and 5, together with a brief discussion of the J/ip polarization. 
Alternative possibilities to produce the J/ip are considered in Section 6. One such mechanism 
is to create a colour-octet cc pair, which then transforms to a colour singlet by rescattering via 
gluon exchange. Another is the associative production of a J/ip and cc pair. The yield from this 
latter possibility is small. However the first mechanism, where two gluons in a symmetric colour- 
octet t-channel state, {gg)ss, belong to two different Pomerons, may dominate at asymptotically 
large energies. The energy and rapidity dependence of J/ip (and T) production is given in 
Section 7. Section 8 contains our conclusions. 



2 The lowest-order amplitude 

We compute the matrix element of the hard subprocess from the diagrams of Fig. |21 in the 
LO collinear approximation. Thus incoming particles, with momenta pi and P2, are taken to 
be on-mass-shell, transversely-polarized gluons. The calculation of the amplitude shown in 
Fig. |2t is similar to the computation of the amplitude for diffractive J/ip photoproduction [H]. 
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Due to the non-relativistic nature of the J /if) wave function, the difference between the quark 
momenta k and k' is very small, that is \k — k'\ <C m. Therefore, following Ref. jjj, we may take 
k = k' = Q/2, and include the integration over the quark loop momentum in the cc — > J/if> 
coupling g, which is normalized to the width of J/ip — > e + e~ decay, see (J3J). 

We expect that the main contribution will come from the region where the rapidity difference 
between the final gluon p and J /if) meson is rather large; that is s = (pi + p^) 2 ^> M 2 . In this 
limit, the amplitude A a , corresponding to the diagram Fig. la, is 

J m A' Nc ^c /> Tr[^(^/2 + m) ^ 2 (-j£/2 - (+ m)^ 2 (j£/2 - p\ + m)] 

where e and e are respectively the polarization vectors of the J / if> and the gluon with momentum 
Pi, and q = p\ — Q is the momentum transferred through the pair of t-channel gluons. N c d abc /8 
is the colour factor, where the indices a, b, c are the colours of the two incoming gluons and the 
final gluon. Bearing in mind possible confinement effects (and to avoid the logarithmic infrared 
singularity as q — > 0) we introduce a cutoff (or effective gluon mass) X g in the denominator of 
(HJ). Since the two t-channel gluons are in a symmetric colour-octet state, the amplitude A has 
positive signature. Therefore it has a small real part, KeA -C ImA 

To obtain the whole amplitude, A, we have to add the contribution of the diagram shown 
in Fig. Et>, which is of the form 

N °**c f,r2 ^5/ 2 Tr[^/2-f+m)^(-^/2 + m)# 2 (^/2-f-^ 1+ m)] 



l mA » = -tS d ahc J dl 2 g(A7ra s f 2 



2tts[(Q/2 -l- Pl ) 2 - m*] [P - X 2 } [(q + I) 2 - Xj] ' 

(5) 

where the minus sign reflects the 'negative' colour charge of the antiquark. We also have to 
account for the contributions where the gluon p\ couples to the antiquark, and not to the quark. 
That is, we must include the graphs with the opposite direction of the 'arrows' in the quark 
loop. Thus the differential cross section is 

d ± = J41 (6) 
dq 2 t 16tts 2 ' K ' 

where the total amplitude A = 2(A a + A b ). 

For transverse (with respect to the p\ il ,p2 i i plane) J /if) inelastic production the trace Tr[...] = 
s 2 m(e ■ e), while Tr[...] = s 2 m(e ■ Q t )/M when the vector e corresponds to a longitudinally 
polarized J/if> meson. Thus, after averaging over the incoming gluon transverse polarizations, 
e L , the ratio of longitudinal (<r L ) and transverse (a T ) cross sections becomes 

da L /dq 2 \Q 2 t \ 



da T /dq 2 t 2M 2 

3 The prompt J/ip yield 



(7) 



In the LO collinear approximation, the cross section of inelastic prompt J/if> production is of 
the form 

da _ rdx 2 , da{s,q 2 ) 
dy^Q 2 ~ J l^ Xl9{Xl)X29{X2) ^qr ' { ) 



where y is the centre-of-mass rapidity of the J /if) meson, x\ = s/ (x 2 S), S is the initial hadron- 
hadron energy squared and the densities of gluons in the incoming hadrons (i = 1, 2). 

For a fixed rapidity of the J/if> meson, the integral over x 2 is equivalent to the integration over 
the mass y/s of the ( J /if) + g) system, where s — {p\ + p 2 ) 2 - 

The hard subprocess cross section a includes the contribution of the diagrams Fig. |2j where 
the single gluon p\ comes from the beam side, plus the 'inverse' contribution, in which the single 
gluon comes from the target. There is no interference between the original (Fig. |2J) and the 
'inverse' amplitudes, due to the different colour structure of single and double gluon exchange 3 . 
Note that single and double gluon exchange correspond, respectively, to antisymmetric and 
symmetric colour octets. 

In the small x region, and at relatively low scales, the gluon distribution behaves like 
xg(x) oc x~ A , where the power A ~ 0.2, while the 'hard' subprocess cross section a does 
not depend on s for s M 2 . Thus the integral over x 2 takes the form J dx 2 /xl +x . The 
main contribution comes from the lowest values of x 2 — M±e~ v /\/~S, which correspond to 
X\ ~ Mj_e y /\/S, where M± = (M 2 + |Q 2 |) 5 . However the essential interval of integration 
available at collider energies ( A In x 2 ~ 1/ A ~ 5) is quite large. This large integration interval 
partly compensates for the small loop (a s ) factor in the g(gg) — > J /if) + g amplitude obtained 
from the diagrams of Fig. EJ 

Moreover for an inelastic process we have to allow for the emission of additional (secondary) 
s-channel gluons from the symmetric octet (gg) (see Fig. EJ). This leads to a power growth of 
the 'hard' subprocess cross section, <r(s), as the function of subenergy yfs. The power behaviour 
of a(s) ~ s 6 is driven by the intercept of the amplitude with the four t-channel gluons (the so- 
called quarteton). Within the leading logarithm approximation (J2 n C n {<^s In s) n ) this intercept 
was evaluated in [Oj. It was shown that numerically the value of 5 is, to within 10%, equal to 
that for the usual BFKL Pomeron (two-gluon) exchange. As a consequence we have to extend 
the dx 2 /x 2 integration over the whole kinematically available rapidity interval Ay. Then the 
prompt J /if) cross section becomes 

- xig(xi)x 2 g(x 2 )Ay da ^^ t - , (9) 



dydQi dqf 

with Xi j2 = (M±/\/S)exp(±y), and Ay = h\{x 2 nSK S/M^). We take the same scale, M±, 
for both gluons; different scales are equivalent to a NLO correction to the LO formula. We 
introduce the factor x max = 0.3 to exclude the contribution from the graphs in which the third 
gluon couples to partons with large Xi j2 > 0.3 in the proton fragmentation regions, that is to 
allow for the fact that parton densities at large x are kinematically suppressed. 4 

The formulae in Eqs. and together with the amplitudes of (|4|5|) . enable the inelastic 

3 The interference would correspond to odderon, instead of Pomeron, exchange in the diagram for the cross 
section. It is known, both experimentally and theoretically, that the odderon-nucleon coupling is small. 

4 The variation of a; max induces a NLL correction; that is, it is equivalent to an a s term (without a ln(l/a;) 
factor) in the BFKL amplitude. Indeed, in the BFKL expansion, the contribution of Fig.^J) to the cross section 
behaves as an a s Ay term as compared to that of Fig.^,. Thus the variation of a; ma x corresponds to an a s term 
without a hi(S/M 2 ) factor. This does not mean, of course, that the variation of x max reproduces the whole 
NLL contribution to the BFKL intercept. 
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Figure 3: A contribution to inclusive prompt J/ip production accompanied by gluon emissions 
from a t-channel gluon-pair in a colour-symmetric octet. Again, the hard subprocess is shown 
by bold particle lines. 
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J /if) cross section to be predicted. In this way, we find that the cross section at the Tevatron 
energy y/S = 1.96 TeV is 5 

a(\y\< 0.6) ~ 2.7 (pQCD estimate) (10) 

in the central rapidity interval \y\ < 0.6 (integrated over the transverse momentum Q t ). Here we 
have taken 6 \ g = 0.8 GeV. This prediction, (JTUJ), is to be compared with the lastest experimental 
measurement [T3] 

a(\y\ < 0.6) = 4.1 ±gj /ib (CDF experiment) (11) 

The uncertainties of the prediction, (|10|) . are as follows. 

(i) The choice of X g . If, instead of 0.8 GeV, we were to take X g equal to 0.5 or 1 GeV then 
u(\y\ < 0.6) becomes 4.0 or 2.0 /ib respectively. 

(ii) The choice of the factorization and renormalization scales. For different scales fi = fiR = 
lip of Mj_/2, Mj_ and 2M±, we obtain cr(\y\ < 0.6) = 2.7, 2.3 and 1.5 /ib, respectively. 

(iii) An unknown K-factor to account for NLO and higher pQCD corrections. 

(iv) The uncertainty in the incoming gluon distribution, which is not well constrained at 
low x and rather low scales. 

(v) The choice of the cut-off x max . The variation of the value of x max plays the role of NLL 
corrections in the BFKL approach. 

Taking all these into account, the expected accuracy of the prediction is about a factor of 
2—3 in either direction or even worse. 

4 Transverse momentum distribution 

Unfortunately we cannot use the amplitudes of Eqs. ()4I5|) directly to calculate the Q t distribu- 
tion of the produced J /if) mesons. In the case of diffractive J /if} photoproduction it is known 
that the interference between the two lowest-order diagrams, Fig. |2k and Fig. |2b, leads to 
a dip in the Q t distribution. Indeed, the differential cross section goes to zero at Q t = M 
[THE)]. However this dip disappears after one includes the leading logarithmic (a s \nS) cor- 
rection The problem is that we cannot resum the analogous corrections in our case, 
since the corresponding 'quarteton' eigenfunctions (coming from 4 gluons in the t-channel) are 
not yet known. 

Therefore we consider a very simple parametrization 

doldQ\ oc g 2 a s (M ± ) 5 \og(x 2 m ^S/Ml)/Ml (12) 

5 We use the LO MRST2001 ^U] gluon distribution at scale /i = M±/2, with the corresponding LO (one 
loop) QCD coupling a s with A^L D = 220 MeV. 

6 The effective gluon mass X g , which occurs in the amplitudes, was first estimated in Ref. where it was 
introduced to describe the photon spectra in J /if> — * jgg decay. A recent evaluation, together with a collection 
of previous determinations, is presented in Table 15 of Ref. ^2). Based on this Table we choose A g = 0.8 GeV, 
with a possible uncertainity covered by the interval X g = 0.5 — 1 GeV. 
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motivated by dimensional counting, which accounts for the dimension of g 2 ~ TM. Again we 
take x max = 0.3. The distribution (fT2j) is normalised by equating its Q\ integral to that of ©• 
In this way the effective gluon mass, A 9 , enters the calculation. The result, shown in Fig. is 
in reasonable agreement with the Tevatron data ■ 

In comparison with the colour octet model, where the J /if) is created in the fragmentation 
of a gluon jet and the expected distribution is da/dQ 2 oc l/Qf , in the case of the Fig. |3] 
subprocess the Q t distribution (JT2J) at large Q t is steeper, da/dQ 2 oc l/Mj_. This is similar 
to the distribution in the colour-singlet model lj. However in contrast with the colour- singlet 
mechanism (Fig. [IJi) here the hadronic transverse energy flow, which compensates the Q t of 
the J/i[) meson, is distributed over a larger rapidity interval in the form of a larger number of 
gluonic minijets (see Fig. EJ). 

According to (J7J) we expect the J/ip mesons to be transversely polarized at small Q t , and 
longitudinally polarized at large Qt, that is Qt 3> M. 7 At present the data are only available in 
the interval of Qt ~ 4 — 20 GeV, where the observed J /if) is approximately unpolarized. The 
parameter a < 0.3 of Ref. [TH] corresponds to a small transverse polarization at lower Q t values. 
However, as Q t increases a changes sign and for Q t > 15 GeV clearly indicates longitudinal 
polarization of the J/if). This is qualitatively consistent with our expectations. In contrast, 
colour-octet models of prompt J /if) production lead to transverse polarization at large Q t [T§j . 

5 Prompt ip f (2S) and T production 

The above formalism can be applied to the production of other quarkonium states by simply 
changing the mass and width of the J p = 1~ heavy QQ resonance. Thus, without any free 
parameters, we can predict the cross section for inclusive prompt quarkonium production. The 
results for if)' and for the upsilon states are compared with the Tevatron data jTHJ in Fig. H] 
and Fig. Irrespectively We use the effective gluon mass X g and the scale that were appropriate 
for the description of the J /if) data. In this way we remove a large part of the uncertainty in 
the prediction for if)'. The good agreement with the if)' data should therefore be regarded as 
support for our perturbative QCD approach. 

The comparison of the prompt T predictions with the Tevatron data, shown in Fig. El is 
complicated, however, since only half of the total T(IS') yield arises from prompt production 
[2~T] . Moreover, the Xb states have a large branching fraction of radiative decays to T(2S); about 
twice as large as those to T(1S). Bearing in mind these complications, and the uncertainties 
in the predictions, the agreement with the data is better than may have been expected. 

The predicted cross sections for prompt J /if), if>'; T(1S,2S,3S) central production at the 
Tevatron energy, = 1.96 TeV, are 

dv/dy\ y =o = 2 - 2 ' °- 6 ^ 40 ' 12 ' 9 nb ' (Tevatron) (13) 

respectively; and correspondingly at the LHC energy, ^/S = 14 TeV, 

do/dy\ y=Q = 8.1, 2.5 fib; 310, 100, 80 nb. (LHC) (14) 

7 At next-to-leading order the polarization of the J/ 1\) may also be affected by the contribution of the longi- 
tudinally polarised incoming gluons (pi in Fig.|2Ji), see for example Ref. [T7]. 



9 




Figure 4: The transverse momentum (Q t ) distributions of inelastic J/ip and ip' production. 
The data are from Refs. US] • The upper and lower data sets for the Q t distribution of the 
J/ip correspond to the total (at = 1.96 TeV) and prompt (at = 1.8 TeV) J/ip yields 
respectively; recall that our QCD prediction is for prompt production only. 
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Figure 5: The Q t distributions of the inelastic production of the T states, compared with 
Tevatron data |2U] . 
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Figure 6: Prompt J/^ (a) inelastic and (b) single diffractive production arising from an inter- 
action in which the two t-channel gluons belong to different Pomerons. One Pomeron is shown 
in bold, the other in fainter print. 



6 Other production mechanisms 
6.1 Uncorrelated gluon-gluon pairs 

Besides the diagrams shown in Fig. |21 where two t-channel gluons in a symmetric colour octet 
state are placed rather close to each other in the impact parameter (bt) plane, there may also be 
a contribution from diagrams like those shown in Fig. El This may be viewed as the production 
of a colour-octet cc pair that subsequently changes colour via rescattering. Here the two t- 
channel gluons in the amplitude belong to two different Pomerons, that is to two different 
parton showers. A convenient way to calculate such a contribution is to use the AGK cutting 
rules |22] , that is to calculate diagram Fig. \§jp and then use the relation 

a FiM = 2 <j Fi ^. (15) 

Strictly speaking, the two gluons in the lower (left or right) parts of Fig. IHb cannot form a 
colour-singlet (Pomeron) state. On the other hand, these two gluons are in a colour symmetric 
state, and therefore we may use the AGK relation (J15)) to simplify the calculation. 

The amplitude corresponding to Fig. IHb is very similar to that of Fig. El The main difference 
is that now the density of the second t-channel gluon is given by an independent gluon distri- 
bution xg(x). That is, we need the probability to find two gluons w gg ~ (xg(x)) 2 . On the other 
hand these two gluons are uniformly distributed over the whole transverse area occupied by a 
proton. Therefore the integration over the transverse momentum Q t is limited by the 'elastic' 
slope B, which may be taken from the slope observed for diffractive J/ip photoproduction; that 
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is B ~ 4.5 GeV 2 |23j. Thus we can consider the lowest-order amplitudes given by Eqs. (|4I5|) 
at Q t = 0, and take the integral / dQ\ = 1/B. 

The factor (A^ c a s /7r) 2 , together with two logarithmic integrations, dl 2 /l 2 , in each amplitude 8 , 
becomes the gluon distribution w gg = (xg(x)) 2 in the cross section. Next we have to account 
for an extra 1/2 in the colour factor, which is cancelled by the factor 2 in the AGK relation 
(|15|). Thus, finally, we obtain 

da Fi ^ Wn^g 2 , , . . , 

— d~ — = 3BM 6 Xl 9{ x V x 29{x 2 ) [xigixi) + x 2 g(x 2 )\ ■ (16) 

Using MRST2001 LO gluons T0|], we find that 9 a Fi ^{\y\ < 0.6) = 2.2 fib, at the Tevatron 
energy y/S = 1.96 TeV. This is smaller than the main contribution discussed in Section 3. 
However at larger energies the 'two Pomeron' cross section will grow faster than the contribution 
of Section 3. Indeed 

a Fi ^ oc x 1 g(x 1 )(x 2 g(x 2 )) 2 oc x^ x xf x oc S 3X/2 (17) 

in the central region of small y, where Xi oc l/y/~S, whereas the cross section © is proportional 
to S x In S only. In particular, at the LHC energy y/S = 14 TeV we expect in the centre of the 
rapidity plateau da Flg ^(y = 0)/dy = 6.7 fib, while the cross section Q is da Flg ^(y — 0)/dy — 

8.1 fib. 

At first sight, the transverse momentum distribution of the J/ip for the contributions of 
Fig. El should be peaked at low Q t , namely Q 2 ~ 1/B. This is true for the contribution of the 
diagram Fig. \§jp, but not for the inelastic process, Fig. The emission of intermediate gluons 
spreads out the J/ip distribution from the inelastic process, so that we expect a spectra 

da Fi ^/dQ 2 oc g 2 a s (M ± ) 3 /Ml, (18) 

analogous to (JEJ). 

We must of course take care of possible double counting. Note that the amplitudes shown 
in Fig. ink and Fig. El are very similar, and the gluon densities given by the global parton 
analyses do not distinguish the gluons coming from one or more parton showers. Therefore, 
to be conservative, in what follows we consider only the contribution of Fig. El A possible 
contribution of Fig. |HK is well within the uncertainties of the lowest-order in a s calculations. 

6.2 Associative (J/ip + cc) production 

Another possibility is to consider the production of a J/ip meson together with a cc pair, as 
shown in Fig. [71 In this case the expected hard cross section at the Tevatron is about a(apc) ~ 2 
nb, leading to 

da(cipc 



dy 



= a Xig(xi) x 2 g(x 2 ) ~ 0.05 fib . (19) 

y=0 



8 The lowest-order gluon distribution at low x is given by the first iteration of the DGLAP equation, that is 
xg{x) = {N c a s /iT) J dl 2 /l 2 

9 We take the scale \i = M/2 both in the gluon distribution and in the QCD coupling a s . 
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Figure 7: The diagram used to compute the associated production of J/ip together with a cc 
pair. The t-channel gluon is needed in order to put the virtual t-channel charm quarks on mass 
shell. 



For inelastic hadroproduction, the contribution coming from this mechanism is not large, of the 
order of 1% . However in e + e~ annihilation, associative production is much more important, see 
for example Ref. [21]. Indeed, in e + e~ annihilation there is no suppression for the production 
of the first cc-pair. On the other hand for time-like annihilation kinematics the third t-channel 
gluon (/ in Fig. |2J) can couple to the nearest parton only 10 . Thus we loose the enhancement 
caused by a large parton multiplicity, that is the factor Ay ~ In (S/M 2 ) in (jUJ). Without this 
factor the contributions of single (Fig. EJ) and associative (Fig. HJ) production are expected to 
be of the same order (in agreement with the measurements of the Belle collaboration for e + e~ 
annihilation [2o]). 

Note that if the inelastic J/ip hadroproduction were to originate from 'colour-octet' dy- 
namics, that is the main yield of J/ip were to come from gluon g* —>■ J/tp + ... fragmentation 
[Hll2E]> then we would expect the same ratio of 'direct' to 'associative' contributions in e + e~ 
annihilation, as that for hadroproduction, contrary to our perturbative QCD predictions. 

6.3 Production via Xc and bb decays 

Finally we have the possibility of non-prompt J/ip production. Experimentally it is observed 
that a fraction of J/tp mesons originate from \ — > J/ip + 7 and b — > J/ip + X decays. We do 
not discuss here the details of these production mechanisms. However in order to predict the 
total J / ip yield we estimate these contributions using the experimental data of Refs. fZT\ 125] 

10 Otherwise we destroy the leading logarithms. 
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and collinear factorization. Symbolically, the relation that we use is of the form 11 

da/dy = xig(xi) a X2g(x2), (20) 

where y = |ln(xi/x2) and the normalisation of o is adjusted to fit the \ and bb data. Since 
the energy dependence of the cross section is driven simply by the x dependence of the gluon 
densities, the experimental data at one fixed energy are sufficient to estimate these non-prompt 
contributions at other energies 12 . 



7 Energy and rapidity dependence 

In comparison with the usual perturbative QCD colour- singlet (Fig. ^i) and colour-evaporation 
mechanisms, where the energy dependence is essentially determined by the product of the gluon 
densities, see (l2Ti|) . our perturbative contribution to the cross section (Fig.HJo) is enhanced by an 
additional logS" factor, Ay in (jHJ). Therefore the energy dependence of these two contributions 
is different, as can be seen in Fig. |HJ To demonstrate the uncertainty arising from the choice 
of the cut-off x max in the additional logarithm Ay in the computation of the new perturbative 
contribution, we show predictions for the two values, x max = 0.1 (lower curve) and 0.5 (upper 
curve). 

Note that the amplitude of the new subprocess is calculated at lowest order in a s . It is 
not unusual to have a next-to-leading order K factor of the order of 2 for the production of 
relatively low-mass states, as in Drell-Yan production for example. On the other hand, we see 
no reason for a strong energy dependence of such a K factor. So, noting the good agreement 
with the Tevatron data, we expect our predictions at the LHC energy to be reliable. 

The energy behaviour of the total inelastic J/ip cross section in the fixed-target/ISR to 
Tevatron energy interval is shown in Fig. El together with the components from the following 
subprocesses: g(gg)ss ~^ J/'i'i 99 "~ > J/^P 9-, Xc J/^ and bb — > J/ijj. We see that the 
cross section is slightly underestimated at the lower energies, where a contribution initiated by 
the qq subprocess, that is by secondary Reggeons, may have some influence. Indeed, it was 
noted in Ref. that the difference between the J/ip production cross sections from pp and 
pp interactions indicates a noticeable qq contribution at the lower energies. 

In Fig. [TUlwe show predictions for the rapidity distributions of J/ip and T(IS') production 
at both the Tevatron (1.96 TeV) and LHC 13 (14 TeV) energies. In Tabled we compare our 

11 The relation also happens to hold for \i production, despite the fact that it cannot be formed by gluon- 
gluon fusion. The dominant mechanism for xi production is the process with antisymmetric colour-octet gg 
exchange. However, since this exchange corresponds to gluon Reggeization, it has the same energy (and x) 
dependence as the exchange of a single gluon. 

12 In Ref. the cross section was given for \ production in the forward hemisphere. Therefore to normalize 
the J/tp cross section, we use l{2"5)) to calculate the integrated cross section. 

13 Note that the General Purpose Detectors ATLAS and CMS can only measure the high transverse momentum 
tail of quarkonium production, the majority of the final-state leptons falling below the trigger thresholds. The 
ALICE detector, on the other hand, is ideally suited to a measurement of both charm and bottom quarkonium 
cross sections © 02] . ALICE can measure electrons and muons down to very low transverse momentum 
(0(1 GeV/c)) in the pseudorapidity ranges \r] e \ < 0.9 and 2.5 < \ri^\ < 4.0 respectively. This gives non-zero 
acceptance for the J/ip down to pt(J/iP) = 0, see Figs. 60 and 63 of Ref. 0. 
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Figure 8: The energy dependence of prompt J /ip production obtained from the colour- singlet 
mechanism, gg — > J ftp g, is shown by the dashed curve. The contribution obtained from our 
subprocess, g(gg)ss J/ip, is shown by continuous curves with x max = 0.1 (lower) and 0.5 
(upper). Also shown are the values for the total J/ip yield measured at RHIC j2H| and the 
Tevatron [T3] . 
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da(pp-»J/y)/dy fib (y=0) 
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Figure 9: The energy dependence of inelastic J /if) production compared with a selection of 
the available data from fixed-target pN to Tevatron pp collider energies pHH |23 Q3j- Besides 
the total (bold curve), we also display the prompt (g(gg)ss — ► J /if) and gg — > (?) and 
non-prompt (% c J/if> and 66 — > J/VO components. 
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pQCD 


CEM 


a j/^ (prompt) (fib) 


80 


33 


<jj/q (total) (fib) 


150 


54 


0"r(iS) (prompt) (/ib) 


2.5 


0.4 



Table 1: The predictions for J/i/j and T(IS') cross sections (in fib) at the LHC For the J/ip, 
the 'total' cross section includes the additional contributions from \ c and 66 decay. The CEM 
predictions are taken from Tables 9 and 10 in Ref. jSj. 

results for the total J/ip and T(IS') cross sections at LHC with the predictions of the CEM 
(taken from Ref. 0). Note that our predicted cross sections are systematically larger because 
of the steeper \/~S dependence caused by the \og(S/M±) factor in Eq. ©. 

8 Conclusions 

We have calculated the prompt hadroproduction of J/ip, ip' and T(1S,2S,3S) states within 
a perturbative QCD framework, without any non-perturbative contributions (such as occur 
explicitly in the NRQCD colour-octet model and implicitly in the CEM). Recall that the original 
colour-singlet LO perturbative contribution (based on the subprocess gg — ► J/ip g) falls well 
short of the data. However, here, we have studied another perturbative QCD contribution, 
which turns out to be dominant. The basic subprocess is g(gg)s s J/ip, see Figs. 12131 or 
Fig. lib- This contribution is enhanced, particularly at high energies, since the additional 
t-channel gluon can couple to a large number of parton spectators. 

The uncertainties of such a computation are listed at the end of Section 3. They are not 
small. However, with our natural choices of scale and of the effective gluon mass, we successfully 
describe the available high-energy RHIC and Tevatron J/ ip data. In addition, without any new 
parameters, we obtain an excellent description of the ip' data, and even a satisfactory description 
of T production. 

There is additional qualitative support for the g(gg)$ s — *■ J/ip mechanism coming from the 
measurement of the J / ip polarization at the Tevatron. This mechanism predicts a longitudinal 
polarization at large Q t , in agreement with the data, whereas the colour-octet model leads to 
a transverse polarization [T9*] . 
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Figure 10: The rapidity distributions of J/ip and T(15 l ) production at Tevatron and LHC ener- 
gies. The continuous and dashed curves correspond to the total and prompt yields respectively 
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